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Abstract: Large-scale coherent magnetic fields are observed in galaxies and clusters, but 
their ultimate origin remains a mystery. We reconsider the prospects for primordial mag- 
netogenesis by a cosmic string network. We show that the magnetic flux produced by long 
strings has been overestimated in the past, and give improved estimates. We also compute 
the fields created by the loop population, and find that it gives the dominant contribution to 
the total magnetic field strength on present-day galactic scales. We present numerical results 
obtained by evolving semi-analytic models of string networks (including both one-scale and 
velocity-dependent one-scale models) in a ACDM cosmology, including the forces and torques 
on loops from Hubble redshifting, dynamical friction, and gravitational wave emission. Our 
predictions include the magnetic field strength as a function of correlation length, as well as 
the volume covered by magnetic fields. We conclude that string networks could account for 
magnetic fields on galactic scales, but only if coupled with an efficient dynamo amplification 
mechanism. 
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1. Introduction 



Large-scale, coherent magnetic fields with strengths in the ;uG range are observed in galax- 
ies and clusters [1, 2, 3, 4]. Models that explain the presence of these fields fall into two 
rough categories: "primordial" mechanisms or "dynamo" mechanisms. The primordial op- 
tion holds that presently observed fields arose from large-scale, homogeneous fields present 
soon after the Big Bang. These primordial fields where diluted by cosmic expansion, then 
slightly amplified by protogalactic collapse to the field strengths presently observed. The 
dynamo option assumes that very weak "seed" magnetic fieds were created by an as-yet- 
unknown magnetogenesis mechanism in the early universe, which were later greatly amplified 
by magnetohydro dynamic (MHD) dynamos operating in spiral galaxies. 

There are many proposals but few truly compelling mechanisms for primordial mag- 
netogenesis. On the one hand, the primordial option provides an explanation of presently 
observed magnetic fields, and the required primordial field strength just evades structure for- 
mation and cosmic microwave background (CMB) constraints [5, 6, 7, 8, 9, 10, 11]. On the 
other hand, it is not clear why homogeneous coherent magnetic fields of the correct strength 
should be present on super-horizon scales. Inflation is one possible cause, but to first order in 
perturbation theory, the vector perturbations required to create magnetic fields decay with 
cosmic expansion. Thus any field generation caused by inflationary perturbations must be 
a second order effect [16]. This initial conditions issue can be circumvented if fields formed 
on sub-horizon scales after the Big Bang, with field strength transferred to the larger scales 
needed to account for galactic magnetic fields [1, 3, 13, 14, 15, 12]. Generally the fields 
produced in this way are very weak, and must be amplified through the action of a galactic 
dynamo or "inverse cascade" arising from turbulent MHD processes [1, 14, 17]. Ultimately 
all of these mechanisms come up against the correlation length problem, which itself arises 
since causal mechanisms can only operate on sub-horizon scales: except for redshifts close to 
decoupling, such mechanisms produce fields on comoving lengths that are too small to ex- 
plain the correlation length of fields observed in galaxies and clusters. For most mechanisms, 
this problem can only be solved by invoking large scale field averaging, inverse cascades, or 
the super-horizon correlations produced by cosmic inflation. The speculative nature of these 
proposals is an indication of how challenging it is to generate fields with the proper length 
scales. In addition, causal mechanisms tend to create fields with a blue power spectrum. If 
the magnetic fields were created before big bang nucleosynthesis (BEN), then conversion of 
stochastic magnetic fields into gravitational waves could lead to violations of BBN bounds 
[18, 19, 20]. 

In this paper we consider primordial magnetogenesis by a cosmic string network. Such 
networks are attractive candidates for magnetogenesis, since long strings typically stretch 
across the cosmological horizon, and so can naturally create coherent effects over large length 

scales. The motion of strings and loops through the primordial plasma also produces the 
vector-type perturbations required for generating magnetic fields. Since the magnetic fields 
are created after BBN, the BBN constraints on stochastic magnetic fields would not apply. 
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For these reasons and others, cosmic string magnetogenesis has been investigated in the past 
by a number of others [21, 22, 23]. These works share a common approach, whereby the 
string network generates vorticity in the primordial plasma, which is subsequently converted 
to a magnetic field through the Harrison-Rees mechanism [24, 25]. The conversion occurs 
because electrons and protons experience differing accelerations due to their Compton scat- 
tering against CMB photons [25]. This differential vorticity in the plasma creates a current, 
which then generates magnetic fields. This process operates so long as the Compton scatter- 
ing of CMB photons from free charged particles is efficient, and so can create magnetic fields 
until decoupling. Conveniently, in this epoch the cosmological horizon is sufficiently large to 
produce magnetic fields on comoving scales which are relevant for the generation of today's 
galactic and cluster magnetic fields. 

We employ a combination of analytic and numerical techniques to estimate the magnetic 
fields produced by the string network. The central challenge is to obtain an accurate estimate 
of the vorticity produced by the network. We argue that the vorticity produced by long, 
horizon-sized strings has been overestimated in the past, and here give improved estimates. 
We also extend our analysis to include the fields generated by the cosmic string loops produced 
in string networks. We find that the loops turn out to give the dominant contribution to the 
total magnetic field produced by the network. In addition to predicting the overall magnetic 
field strength, we use properties of the string and loop networks to compute the spectrum 
of field strength as a function of correlation length, as well as the fractional horizon volume 
coverage of regions in which fields are generated. We obtain these predictions from a computer 
code which implements two semi-analytic string network models, the one-scale model (OSM) 
[26] and the velocity-dependent one-scale model (VOS) [27, 28]. The code also incorporates 
the dynamics of individual loops, including the forces and torques on the loops from Hubble 
expansion, dynamical friction, and the emission of gravitational waves. We make minimal 
assumptions about MHD processes, and do not invoke any large scale averaging, inverse 
cascades, or turbulence. 

We find that for reasonable choices of string network parameters, the network can create 
sufficiently strong seed fields on galactic scales, provided the subsequent dynamo amplifica- 
tion is very efficient. Cosmic string magnetogenesis has the satisfying property that it makes 
a number of other predictions. We find that adequate magnetic fields can be generated by 
cosmic strings with tension Gji/c? > 10~^, a value which is still allowed by cosmological 
observations [29, 30, 31, 32, 33], but may be constrained further by an array of other cos- 
mological observations set to be made within the next few years [34, 35, 36, 37, 38, 39, 40]. 
Various loop network models give similar predictions for the magnetic field strength as a 
function of correlation length, and this scale-dependence might be observationally testable. 
Optimistically, the discovery of cosmic strings with relevant tensions would greatly bolster 
the viability of this mechanism. On the other hand, if galactic dynamos turn out to be less 
efficient than we assume, or if no strings with a high enough tension are found, then our 
mechanism is ruled out. 

We organise this paper as follows: in Section 2 we describe our estimates of the vorticity 
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generated by long string encounters, and discuss the differences between our long string esti- 
mates and those in the literature. We also give estimates for the vorticity generated by loops, 
leaving a detailed discussion of the loop dynamics for Appendix A. Armed with expressions 
for the vorticity generated by individual strings and loops, we give a brief description of the 
network models that describe the long string and loop populations in Section 3. Wc then use 
this information to present analytic and numerical estimates of the magnetic fields generated 
by the string network in Section 4. We give the magnetic field strength as a function of 
correlation length, describe the volume fraction covered by the loops, and discuss how our 
results depend on the various string and loop network parameters. Some details of the code 
we developed to obtain these results is given in Appendix B, and wc tabulate our parameter 
choices in Table 2. Our summary and conclusion appear in Section 5. Throughout this work, 
we set c = 1 but keep other natural constants (such as G) explicit. 

2. Generating Vortices 

Long strings generate vorticity by dragging the plasma as they pass through it, and one 
computes the drag force by calculating the impulse given to test particles by the passing 
string. This problem is considered in detail in refs. [41, 42, 43, 22, 23], from which we 
draw some essential facts, summarised below. The simplest situation to study is that of 
an infinitely long, straight cosmic string with tension T, fundamental mass per unit length 
Ho, and effective mass per unit length fi. If extended along the z-axis and moving in the 
x-dircction with velocity v^, it imparts a velocity Vy to a test particle located far from the 
xz-plane, directed toward the x^-planc and given by 

Vy = h 4TrGiJ,oVsjs ■ (2.1) 

where 7^ = (1 — v^)'^^"^, and A is defined as 

X = fi-T. (2.2) 

The effective mass per unit length, /x, and tension, T, are obtained by averaging over the 
small scale fluctuations, or wiggles, on the string. They are related to the fundamental 
"bare" tension /xo by [41, 44] 

fxT = i4. (2.3) 

Typically /j, « 1.9/io in the radiation era and /j, ss 1.5/Lto in the matter era [45]. Equation 
(2.1) is valid when we can work in the weak-field limit of GR, when the test particle is far 
enough from the string that the wiggles can be effectively averaged out, and when the string 
is moving sufficiently rapidly that the test particle does not move substantially closer to the 
string during the encounter. 

Only the first term in (2.1), which dominates for small string drift velocities Vs, is im- 
portant for this work. This term is due to the Newtonian gravitational force of the string: in 
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linearised GR the Newtonian potential is sourced by both energy density and pressure, so the 
string creates the same acceleration field as would an infinite rod with linear mass density A 
in purely Newtonian gravity. The second term in (2.1) arises from the conical deficit angle 
SttG/xq which the string introduces in the surrounding spacetime. Simulations indicate that 
wiggles move very quickly, with an RMS velocity of vrms > 0.6, and that the velocity of 
strings averaged over a correlation length is quite small, Vg = vrms ~ 0.15 [41, 46]. So typi- 
cally the first term dominates, and since 7^ f« 1, we are safe in using the non-relativistic limit. 
Therefore we are justified in treating the physics of string encounters with Newtonian gravity. 
Unlike the conical-deficit approach to long strings, this generalises to loops in a simple and 
convenient way. It is also a conservative approach, for including the deficit angle contribution 
leads to stronger plasma flows in the string wake, enhancing the vorticity so produced. 

Using the deflection formula (2.1) we re-examine the calculation of vorticity created by 
long strings in Section 2.1. We then apply the same method to estimate the vorticity generated 
by string loops in section Section 2.2. We show that producing vorticity by long strings is less 
attractive than previously thought, but that string loops provide a well-motivated alternative 
mechanism. 

2.1 Straight Strings 

Long strings generate vorticity by dragging the plasma behind them as they pass through it, 
or equivalently by producing a velocity component Vx on test particles. Treated as a first-order 
perturbation, the x-component of the force on a test particle vanishes when averaged over 
the encounter with a string. At second order one includes the motion of the test particle in 
the y-direction during the encounter, which brings the particle slightly closer to the xz-plane 
during the departure of the string than it was during the approach, resulting in a net force 
in the x-direction. The final velocity can be computed [21] by going to the string rest frame, 
where the test particle moves with velocity Vg ■ Since the Newtonian force is conservative, the 
magnitude of the test particle's velocity is unchanged by the encounter, so 

vl=vl + {vs-vxf, (2.4) 

and therefore 

vx = ^^+0{vl). (2.5) 

To generate magnetic fields, the passing string must generate vorticity in the plasma. After 
the encounter, matter moves toward the string's trajectory with velocity Vy, but this flow 
is not rotational. The plasma's total momentum vanishes after the encounter, due to the 
symmetry of the flow with respect to the a;^;-plane. At second order, the flow created by 
the drag velocity Vx has net momentum in the direction of the string movement, but is not 
rotational either. 

A rotational flow can be created by vortices that build up due to turbulence in the 
string wake, as first proposed by Vachaspati and Vilenkin [41, 42] (see also [47, 21]). The 
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Figure 1: Two straight strings with effective Newtonian mass density A = /i — T cause wakes in the 
surrounding plasma via gravitational interaction. After the encounter, the magnitude of the dragging 
component of the plasma flow velocity is approximately Vy ^ GX/vg and ~ Vy/(vs). The resulting 
plasma flow carries a rotational component of velocities up to Vrot Vx over interstring distances, so 
that an angular frequency of order Cbpi ~ X^G"^ /{RgV^) results, see (2.9). 



authors of [21] argue that a two-string encounter creates a vortex whose size is comparable 
to the interstring distance, with rotational velocities comparable to the infall velocity of the 
plasma Vrot ~ Vy. Dimopoulos and Davis [22, 23] later investigated the two-string encounter 
assuming the same plasma flow. The argument they employ runs essentially as follows: 
plasma of density /o in a region of volume V ~ R'^VgT, with post-encounter net momentum 
Ap ~ Rgpvx, leads to a force on the string of F ~ Ap/T R^pVy. Applied over a distance 
Rs, the string does work Wg ~ pR^Vy. If we were to assume that the overall flow in the region 
V is rotational and were to set Erot ~ P^s^'^ot equal to Wg, this would imply Vrot ~ 
The keystone of these proposals is the assumption that the vortex size is comparable to the 
interstring distance, and that it carries the majority of the total kinetic energy imparted to 
the plasma. As we have shown, such a vortex is not present immediately after the encounter. 
Indeed, the change in angular momentum of the strings due to the dragging of the plasma is 
roughly AJg r-^ RsAp ~ R'^pv^- Conservation of angular momentum implies that the plasma 
may have, at most, a rotational component in the volume ~ R^ with angular momentum 
Jpiasma ~ R'tp'^rot ~ AJ^. Therefore the net rotational velocity is closer to Vrot ~ Vx, which 
is much smaller than Vy. This is illustrated schematically in Figure 1. 

These estimates can be improved as follows. The wake behind a single string, created 
at ^ teq, has a length l^j, width Ww and thickness dw, given by the scaling relations 
[48, 49, 45, 50] 

i-w ~ tp — , Ww ~ Vgtp — , dw ~ Vytp — , (2.6) 
z z \ z ) 



^An additional factor of 2 in [23] stems from taking Vrot ~ u = 2uy, where u is the relative velocity of 
particles on opposite sides of the string. 
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which are vahd for z > zpVy/vs- Turbulent eddies arise within the wake shock [47], and 
could potentially lead to large rotational velocities of order Vy. However, the characteristic 
size associated with matter chunks due to fragmentation of the wake is comparable with the 
thickness of the wake behind an individual string [52, 51, 45]. Wc expect turbulent, gravi- 
tationally supported vortices at this length scale. Comparing this thickness to the interstring 
distance Rs{t) ~ P^Vgt yields 

where we use Vy from (2.1), P is the intercommutation probability ^, 1/2 < /3 < 1 [46], and 
tp/t = {ap/af^"^ ~ {zlzpfl'^ . Considering mildly relativistic strings [vg ^ 0.1 as an order of 
magnitude) and the largest possible G\ ~ 10""^, we are left with 

Rr^\~^- 

Below we show that the vorticies relevant for magnetogenesis are created toward the end of the 
radiation era, so the redshifting factor is of order unity and vortices due to turbulence in the 
string wake are much smaller than the interstring distance. The seed fields necessary to initiate 
plausible galactic dynamos should be coherent over distances of at least ^seed ~ 5 — 50 pc at 
decoupling, but (i^ ~ Ipc for turbulent eddies created around t^ec- We discuss the relevant 
length scales in more detail in Section 4. 

Since the magnetic field is directly proportional to the angular velocity of the plasma, 
our new estimate greatly reduces the expected strength of magnetic fields produced by the 
long strings. Our arguments indicate that the drag velocity f;^. is the relevant velocity for 
magnetic fields that arc coherent over interstring distances, not the infall velocity Vy. This 
gives a plasma vorticity 



^pl^ ^ ^ r,r, ~ on 3 • strmgs, this paper) (2. 



This is in contrast with [22, 23] which obtains the estimate 

[23] Vy lirXG 

vl ~R ^ — p~ ■ (long strmgs, previous) (2-10) 



We find that while rotational velocities ~ Vy are possible within the string wake, by (2.8) 
their correlation length is much smaller than the interstring distance Rg. 

To summarize: vortices due to turbulence in the string wake are confined to small scales, 
much smaller than the interstring distance Rg- Even though these vortices may have large 
rotational velocities, their small size makes them far less appealing for magnetogenesis. 
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Figure 2: a) A rotating loop with angular velocity uj and drift velocity Vt attracts the surrounding 
plasma with a net velocity of order Vy ~ GX/vt, which in turn causes a vortex with rotational velocity 
of order Vx ^ Vy/vr over the size of the loop. The resulting angular velocity of the plasma is of order 
i^pi ^ X'^G'^ /{£VfVr), see eqn. (2.25). b) Shematic of a rotating loop at rL(o', i) and a test mass at r^. 
Here (j) = ujt, 9 = ira. 

2.2 A Rotating Loop 

We now consider the generation of rotational velocity flows in the plasma by a loop instead of 
a long string. We make the assumption that a loop in a typical state of motion has a nonzero 
angular velocity lo. The drag force from the loop then transfers angular momentum to the 
plasma in a straightforward fashion. The correlation length of the rotational flow, and thus 
of the magnetic field, is set by the size of the loop. In this section we estimate the resulting 
vorticity. 

As in the long string case, our first task is to compute the gravitational impulse exerted 
on a test particle by a passing rotating loop. We consider a rigid, circular loop of radius 
R, length (. = 2itR, and linear mass density A, with its angular velocity cu and translational 
velocity vt oriented along the z-axis. This situation is illustrated in Figure 2. Analytic 
solutions to the loop equations of motion are known [54], but in general the loop dynamics 
are quite complex, and certainly we do not expect cosmic string loops to act precisely as 
the rigid loops we study here. Nonetheless, we assume that the relevant physics operating 
on the largest loop length scales is effectively captured by the idealization of a rigid loop, 
and its "coarse-grained" properties such as velocity, angular momentum, mass and size. The 
parameters i, Vr and vt are all functions of time, thanks to the dynamical forces acting on the 
loop and its emission of gravitational radiation. We describe the equations that govern these 

« 1 for cosmic strings, 10~^ < < 1 for D-strings and 10"'^ < f < 1 for F-strings [53]. 
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quantities in Appendix A, and find that I, Vr and vt change very httle over the timescales 
associated with test particle encounters. While we take all of these dynamical forces into 
account when studying the long-term evolution of the loop population, for the purpose of 
estimating the drag on the plasma, we treat £, Vr and vt as constants. 

Switching to the loop's rest frame, we consider a particle at r^^ = Ry. At time t, we 
parameterise points on the loop by 



/ sin-TTcr cos U!t\ 



rL = R 

where a ranges over —1 ... 1. We take the ratio 



siuTTO" sina;f 

COSTTO" 



(2.11) 



n = ^ (2.12) 

to be larger than, but close to, unity. We further define the displacement d = r/, — with 
magnitude 

d = R(^l-^smTTasmLot+'^y . (2.13) 

Recalling that we have defined x as the direction parallel to the string's motion, we know 
that the acceleration component vanishes to first order if we average over a full rotation. 
The net acceleration is in the y-direction, toward the loop, and given by 



ttGXR I drj da^=Ci7r-^, (2.14) 



where r = t/T, with T the loop rotation period T = 27r/a; and 

Ci = R^ J\t J\a^, (2.15) 

which is of order one. (For example, one finds Ci ~ —1/2 for R = 2R). The net velocity 
toward the loop after one rotation is then 

Vy — \ay\ = \Ci\2tt^ ^TT^ . (2.16) 

Thanks to its drift velocity wt, the loop undergoes roughly 4Rij/(2Trvt) = 2v,-/{TTVt) rotations 
before it moves away from the test particle, so the total velocity acquired by the test particle 
during the encounter is 

vy-^——, (2.17) 
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which is similar to the straight string case in (2.1). Only the translational velocity vt enters 
this expression, since the longer a particle experiences the gravitational attraction toward the 
loop, the faster they approach each other in the end. 

The drifting loop drags the plasma behind it - just as a straight string docs - but unlike 
the string encounter, the flow has a rotational component. The drag velocity in the a;-direction 
is again of order Vx ~ Vy/vt (see below), resulting in a dynamical friction force on the loop of 
F ~ R^pVy ~ R^pG^X'^/vf. This has two effects on the loop. First, it feels a net force due to 
dynamical friction of [55, 56] 

^t = ~lne-l (2.18) 

6* 



where 



t ' max 

is the minimum scattering angle, Vmax = J vtdt o:^ Svt, and 

Second, because the loop rotates, the drag force generates a vortical flow. As in the long string 
case, this is a second-order effect. The acceleration in the x-direction due to an infinitesimal 
element of the loop dcr is 

ttGXR'^ . , 
dax = — ^3 — sm TTcr cos a;t dcr (2-21) 

Substituting the first-order trajectory of the test particle, given by 

ym{t) = R + Vy(t--) (2.22) 
\ (jj / 

with ^^(vr/a;) = R, into (2.21), expanding in terms of e = Vy{t — T/2)/R <C 1, computing 
the drag velocity by integrating the first order term over a single period (the zeroth order 
contribution vanishes due to symmetry) , and replacing GX in terms of Vy results in 

2 2 

~ !i|C2|^-^, (2.23) 

where 



-1 Jo 



C2 = 127r / da / dr sin vra cos 27rT ( ~ 2 



TZ — sin vra sin 27rT 



(1 + 7^2 - 27^ sin ira sin 27rT)5/2 



(2.24) 
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Numerical integration gives C2 ~ —0.14 for TZ = R/R = 2, so we take C2 = —1/7. The drift 
velocity Vt enters (2.23) via Vy lirGX/vt, so as a result of the drag force, the plasma is 
stirred up with angular velocity 

where we take the vortex size to be given by the loop length i = 2itR. For comparison, the 
vorticity due to straight string encounter, derived in Section 2.1, is 

where Rs denotes the inter string distance Rg ~ P^Vst. 

We show in Appendix A that Vr{t) and vt{t) are comparable to the average velocity of 
straight strings, Vg- Therefore the ratio of vorticities is essentially controlled by the ratio £/Rs 
of loop length to long string separation. Since the loop must linger in each region of space 
long enough to establish a vortex, vortices are only created when Vr > Vf. 

In this section we focus entirely on the generation of vorticity, since it is a precondition 
for the Harrison-Rees mechanism to create magnetic flux. Since these vortices are real astro- 
physical objects, they are subject to many physical processes which we have not considered 
here, and our comparisons of the length scales and vortex velocities should be viewed in this 
light. For instance, though there is a small vortical component to the plasma flow in the 
region between two widely separated cosmic string wakes, this vast distance will encompass 
many local over- and under-dense regions, complicating the physics. On the other hand, the 
length scales perturbed by string loops and single string vortices are somewhat smaller, and 
should thus be less subject to the vagaries of small plasma variations. Furthermore, in each 
case the strings are doing more than generating vortical motions and magnetic fields: they 
are accreting matter themselves. This adds to local overdensities, and makes the regions over- 
swept by strings somewhat more likely later to collapse and develop structure. We should also 
have string loops attracted to relatively overdense areas. We may thus expect regions that 
have been affected by string magnetogenesis to be, preferentially, those regions which later 
form structure. Though small, we believe this phenomenon will help to increase the effective 
coverage of string-sourced magnetic fields, since even if they fail to cover the whole universe, 
the parts that they do cover will likely be the parts that will eventually host galaxies. 

3. String Network Models 

To obtain a prediction for magnetogenesis, we must combine our results for the magnetic fields 
from a single loop or string, obtained in Section 2, with a model of the network population. 
We use the model to provide the loop length spectrum, defined by 

iV(At) = ^^ (3.1) 
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where ^ is the loop length and Nioops the number of loops per comoving Hubble volume. To 
explore the dependence of our results on network model assumptions, we use two popular 
semi-analytic loop network models: the one-scale model (OSM) and the velocity-dependent 
one-scale model (VOS). Both models obtain equations for the evolution of the long string 
energy density, and then apply energy conservation to determine the energy fed into the loop 
population in the form of newly formed loops. 

3.1 The one-scale model (OSM) 

According to the OSM [26, 43, 45] the absolute number of loops Nioops in a physical volume 
V{t) obeys the equation 

dNioops ^ V{t) C 

dLnit) Lnity a ^ 

where C is constant during matter- or radiation-dominated epochs, a is the size of newly 
created loops as a fraction of LH{t), and 

LH{t)=a{t) / — - (3.3) 

is the particle horizon measured in physical units. A loop formed at time is taken to have 
initial length (.{tp) = aLffitp), and then subsequently loses energy through the emission of 
gravitational radiation. The length at time t of a loop formed at tp > t is given by 

£{t, tp) = frOiLn {tp) - TeGfMit - tp) (3.4) 

where fr represents the energy loss from the redshifting of the loop peculiar velocity imme- 
diately after formation, and is a dimensionless parameter controlling the efficiency with 
which the loop emits gravitational radiation. 

Using the fundamental OSM equations (3.2,3.4) we now compute N{£,t). To integrate 
(3.2) we take V{t) = {a{t)RY, so (3.1) measures the absolute number of loops per (cubical) 
comoving volume R^. As we are primarily concerned with the radiation dominated epoch, 
we take Lnit) = 2t = 1/H{t). To calculate N{£,t) we first rewrite (3.1) as 

dLnitp) dine ■ ^'^■^> 

By (3.4), a loop currently of length £ at time t must have formed at a time tp when the 
particle horizon Lnitp) was 

which yields the second factor in (3.5) by differentiation. To obtain the first factor, we pick a 
fiducial time to during radiation domination, set a{to) = 1, and fix i? = Lj^(io)- This yields 
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the number of loops per logarithmic interval in ^, at time t, in a comoving volume equal to 
one Hubble volume at a time to, as 

N{i, t} = — [afr + ^— (3.7) 

where we have taken the value C ~ 30 during radiation domination, as in [26]. Despite the 
fact that new loops are being created dX £ = aLn, the spectrum possesses a peak at the 
characteristic length 

ipeak{t) = ^^LH{t)^50-Gfi-ct (3.8) 

taking ^ 50 [26, 79, 80]. The spectrum falls off for i > Ipeak since the large loops are 
created during times of slower loop production, and falls off for £ < £peak since the small loops 
are near the end of their lives and are evaporating rapidly. Similar expressions can be derived 
in the matter era. 

3.2 The velocity-dependent one-scale model (VOS) 

The VOS model of refs. [27, 28] characterises the string population by a length scale L, a 
velocity v, and a string number density n. Commonly, one makes the approximation n = 
but for our code we maintain these as distinct parameters for flexibility. Like the OSM model, 
the VOS model has several dimensionless parameters, which wc call ci, C2, and C3. These are 
fixed by matching to numerical simulations. Taking from [27] the scaling values of HL and 
V, we find ci = 0.21 (0.2475), C2 = 0.18 (0.3675) in the radiation (matter) eras. We allow 
free evolution of the network in the dark energy era using matter era parameters, as this 
epoch does not effect magnetogenesis. The third parameter, C3 = 0.28, fixes the scaling 
string number density. The physical origin of these equations and their parameter inputs is 
discussed extensively in ref. [27]; matching between the equations presented here and those 
of ref. [27] is explained in ref. [28]. The VOS model assumes that the length scale evolves 
according to 

^ = HL + civ (3.9) 



where the loop parameter ci < 1 is dimensionless. The velocity v obeys 

dt; _ -2Hv + C2/L 

and the comoving number density of long strings N = a^n is governed by 

dN C2NV csN'^Lv 



dt L a 



2 



(3.10) 



(3.11) 



To translate the output of these equations into what wc want - namely, loop creation rates - 
requires a further step. Taking poo = 'T'A') we can just plug the rate of change of cosmic string 
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energy density into the loop creation rate formula from ref. [26], 




)) 



(3.12) 



where in this formula Lh is again the horizon size, V is the horizon volume, and a is, again, 
the size of the loops after creation as a fraction of Lh- 

4. Magnetogenesis and Observational Bounds 

We now gather our preliminary results to determine whether the string network can account 
for the magnetic fields observed today. In Section 4.1 we follow the evolution of the fields from 
decoupling, through galaxy formation, to the operation of a dynamo mechanism. In Section 
4.2 we give analytic estimates for the resulting galactic fields, followed by our numerical results 
in Section 4.3. There are great uncertainties affecting almost every stage of the evolution of 
magnetic fields - especially the dynamo amplification efficiency - but we show that with an 
efficient dynamo it is possible to account for presently observed galactic magnetic fields. We 
summarise our parameter choices in Table 2, found on page 37. 

4.1 Evolution since matter-radiation equality 

The vorticity created by the cosmic string network sources magnetic fields through the 
Harrison-Rees mechanism [24, 25] . Harrison [24] was the first to consider a model whereby 
a seed field is generated by a vortex in the ionized plasma. In Harrison's model, a circular 
current develops when there is a differential rotational flow between ions and the electron- 
photon gas. Due to Thompson scattering, non-relativistic electrons are tightly coupled to the 
radiation bath before tdec- Ions are not tightly coupled, so expansion damps their angular 
momentum more efficiently. Because the electrons outpace the ions, a net current is set up 
which generates a magnetic field. Rees [25] pointed out that Harrison's mechanism creates 
vortices that decay too quickly in the radiation era to produce useful fields, and suggested a 
modification. Rees' mechanism relies on vortical fiows in the early matter era, when vortical 
motion is less destabilized by expansion. Since the plasma is still optically thick before recom- 
bination, the angular velocity of the electron fluid is damped through Compton drag caused 
by the microwave background; meanwhile the relatively heavy ions are less affected. Again, 
a current is established, though in a direction opposite that of Harrison's original proposal. 
The strength of a magnetic held produced in this way is [21] 



where B is in Gauss and copi in s . This mechanism can only produce magnetic fields 
when the universe is ionised and Compton scattering is efficient, which ceases to be true 
after decoupling. Though it is also true that the universe reionised at low redshift, the 
Harrison-Rees mechanism no longer works in the late universe since the radiation density is 
low, rendering Compton scattering ineffective. 



B = LOpi 10 CJpi , 



(4.1) 
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Once magnetic fields are produced at a redshift zp, their (proper) correlation lengtli ^ 
grows with the expansion of the universe as 

and their field strength evolves according to 

as a consequence of flux conservation. Once galaxy formation begins, the evolution of mag- 
netic fields becomes far more complicated. As a protogalactic cloud becomes non-linear and 
begins to collapse, the correlation length decreases but the field strength is amplified. While 
there may be some amplification of the field during collapse [57, 58], we assume that no dy- 
namo is active at this stage and the field strength is primarily governed by flux conservation 
[59, 60]. The net predynamo amplification factor in a spiral galaxy is approximately [15]^. 

B 

8 X 10^ (4.4) 



^9f 

where Bgf = B{zgf). Once the protogalactic cloud collapses, we assume that the field is 
amplified exponentially by a dynamo mechanism [15], such as the aw-dynamo ^. This dynamo 
begins operation when stellar winds and explosions generate interstellar turbulence, which 
transforms into cyclonic motions through the Coriolis forces associated with galactic rotation 
[61, 62]. The magnetic field surrounding the galaxy has two modes, a toroidal and a poloidal 
component. The dynamo converts the poloidal to toroidal flux by differential rotations of 
the galactic disk (the cj-effect) and the toroidal to poloidal through the cyclonic motions (the 
a-effect) [63]. The combined effects can amplify the magnetic field strength by many orders 
of magnitude [64, 65, 66, 15]. 

We parameterise the dynamo by an efficiency factor Tdy, such that V^^ is the field strength 
e-folding time, and the field Bq measured today is related to the initial field Bi by 

ln^ = rd,(t/-i,) (4.5) 

where U > tgf indicates the onset of dynamo amplification, tf < to is the time at which the 
fields reach the observed value, and tgf is the time at which the protogalactic cloud collapses. 
We take ti = tgf and tf = tQ to arrive at the most optimistic lower bound of the seed field. 
The resulting amplification is weakly sensitive to the choice of tgf , and here we take tgf = 475 
Myr, which corresponds to Zgf = 10. 

Although the value of Fdy is crucial to estimate the necessary seed field, its value is a 
matter of considerable debate. In the literature, one finds many values for T^y scattered in 



^This expression includes: formation of a halo, gain of angular momentum through tidal interactions with 

neighbouring galaxies, and disc formation, sec [15] for a review. 

''The aw-dynamo assumed here applies to galactic magnetic fields growing in discs only. 
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galaxy formation 






0.2 Gyr 


0.3 Gyr 


0.5 Gyr 


1.0 Gyr 


6 


IGyr 


3.8 X 10^^ 


2.4 X 10^8 


1.1 X 10^1 


3.3 X 10^ 


10 


475 Myr 


(*) 5.2 X 1028 


1.4 X 10^9 


3.1 X 10^1 


5.5 X 10^ 



Table 1: The amplification factor Bo/Bi, tabulated with a variety of assumptions regarding the time 
tgf and redshift Zgf of galaxy formation and the efficiency Tdy of the galactic dynamo. The value we 
take to get the most optimistic lower bound on Bgeed in (4.6) is marked with a 

the range 0.2 Gyr < T^^ < 0.8 Gyr [68, 64, 67]. Recently, some have contended that even 
larger values, T^J^ > 1.1 — 1.4 Gyr, are more likely [15]. Still, taking this uncertainty as an 
opportunity for optimism, we consider very efficient dynamos with V^^ = 0.2 Gyr, which 
amplify by a factor 5.2 x 10^^ (the amplification factors for various choices of T^y and tgf are 
given in Table 1; note that in the numerical results that follow, §4.3, T'^^ = 0.3 Gyr will be 
adequate). To obtain the present field of Bq = 10^^ Gauss under these assumptions, the field 
Bseed that must be present at decoupling with Zdec = 1089 is 

Bseed ~ 2 X 10-35 G . (4.6) 

This is a very optimistic lower bound for the seed field at tdec- For seed fields between this 
limit and lO^^*^ G, only the most efficient dynamos might work, though the existence of such 
dynamos in nature is controversial [15]. 

In addition to constraints on its strength, the seed field must possess a sufficiently large 
correlation length. The correlation length after protogalactic collapse ^gf must satisfy ^gj > 
100 pc for the dynamo to commence [60]. Using a simple spherical collapse model for galaxy 
formation to estimate the comoving correlation length Xcorr before galaxy collapse leads to 
[60] 

Xcorr > V Xgal = 0.95 Tj (O^/l')-^^^^?!^^ [Mpc] (4.7) 

where Mu = M/IO^^Mq, M is the mass of the galaxy, Xgai is the comoving length of the 
galaxy at formation, and rj is the fraction of a galaxy over which the magnetic field has to be 
correlated. Taking M12 ~ 0.1 and rj pa 1/150 (corresponding to ^gf 100 pc) yields 

Xcorr ~ 5.8 kpc. (4-8) 

Gonsequently, the seed fields must have a physical correlation length ^seed at decoupling of 

6eed=T5^>5-4pC. (4.9) 
1 + Zdec 

which compares favorably with the particle horizon of f« 200 kpc at z f« 1000. Larger seed 
field correlation lengths are even better: the 5.4 pc minimum quoted here should cover only 
about 10% of the protogalactic cloud, which is only marginally adequate. A seed field with 
a longer length scale - say 50 pc - would comfortably suffuse the whole protogalactic cloud 
with a single coherent field. 
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4.2 Analytic Estimates 
4.2.1 Near decoupling 

Straight string encounters near decoupling produce vorticity given by LOpi from (2.26), which 
by the Harrison-Rees mechanism creates a seed field of 

Bs ~ W-^{27r)'it^^ll^ < 1.6P-^ X 10-26 G, (4.10) 

where we use G(/Lt — T) < 10~^, Vg > 0.1 and an interstring distance of Rg = P^Vgtdec- These 
fields have a correlation length at decoupling of 

~12P^kpc. (4.11) 

The field strength in (4.10) is several orders of magnitude smaller than the corresponding one 
in [23] , since our estimates of the vorticity on large scales generated by a two-string encounter 
is much lower than that in [23]. Cosmic strings with ~ 1 and the largest possible string 
tension produce fields that are just strong enough to seed the most efficient dynamos. This 
improves somewhat for F and D-strings, which can have a lower P^. In either case the 
coherence length is larger than the minimal one in (4.9). 

The vorticity from rotating loops with Upi from (2.25) results in seed fields of 

Taking G{ji — T)< 10"^, vt > 0.1, Vr ~ 0.4 and defining a new parameter a such that the 
loop length ^ at formation is 



tp = tdec we obtain 



B,<'-^^G. (4.14) 

a 



We use Vt > 0.1 since even a large initial velocity of vt ^ vrms decreases due to redshifting 
in the matter era, before it speeds up again due to the rocket effect (see Appendix A. 2). 
Because redshifting is absent for the rotational movement, the rotational velocity decreases 
only very slowly due to the emission of gravitational waves, which is also counterbalanced by 
loop shrinking to some extent (see Appendix A.l); hence we use use n,. ~ 0.4 (Appendix A. 3). 
Since a < P^Vg, we can achieve a larger field strength than for a straight string encounter. 
Consequently, less efficient dynamos work, but not all galaxies are so lucky as to have had a 
loop sweeping over them in the past. The coherence length 

^i~117akpc. (4.15) 
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is large enough to seed the dynamo for de > 10~^, so that both the largest loops and many 
of the smaller ones contribute to magnetogenesis. The smaller the loops, the stronger the 
resulting seed field. 

The analytic estimates suggest similar contributions in magnetic field flux from long 
strings and loops. In fact, the numerical estimates presented in Section 4.3, which include 
more details of the loop dynamics and population characteristics, show that loops produce 
much stronger magnetic field fluxes than long strings. Partly this is because the loop length 
spectrum peaks at lengths much smaller than aLn, which effectively lowers the value of a 
and greatly increases the fields they create. In addition, redshifting slows the loops, which 
then create stronger magnetic fields. 

4.2.2 Fields from before Matter domination 

Extending magnetogenesis into the era before matter domination is a tricky buisness. As 
one pushes to higher redshift, the horizon becomes ever smaller, and sufficient correlation 
lengths even harder to achieve. Also, having more expansion time for vortices and magnetic 
fields to dilute pushes the predicted field strengths from these earlier times to even lower 
and less plausible values. In addition, many dissipative processes are present in this era, 
causing both B and ^ to decrease further; as a consequence, it is unlikely that any primordial 
magnetic fields within the horizon survive with reasonable strength and coherence length [70] . 
To complicate matters even more, the full set of MHD equations should be used deep in the 
radiation era [69], instead of the simple scaling laws. We conclude that magnetic fields with 
the largest fluxes and correlation lengths are the ones created around decoupling. 

4.3 Numerical Results 

The mechanism for magnetogenesis presented here involves too much physics to be treated 
efficiently using only analytic methods. But each of the relevant physical effects is reasonably 
well studied, so we include them using a computer code. The code is discussed in more detail 
in Appendix B. Briefly, the advantages of using the code include: 

• Loops are dynamical objects, and a code can track their evolution over time. 

• The fraction of the universe's volume that is given a seed field is very important to 
know, yet very difficult to estimate analytically. In a code, this can easily be computed. 

• String parameters - like the bare tension Gfj^o, and a, the loop size - enter the equations 
in a variety of places, making it hard to guess how magnetogenesis strength, length 
scales, and volume coverage will be affected by each. Our code allows us simply to try 
a variety of values and compare the outcomes to develop some intuition about their 
effects. 

• The relative sizes of loop and straight string contributions to magnetogenesis are easy 
to calculate and compare on equal footing, without making the possibly prejudicial 
approximations necessary for analytical estimation. 
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The model parameters we wish to vary are: 

1. GjJsQ, the string's bare tension; 

2. our model for string network evolution; 

3. whether or not string loops are allowed to undergo dynamics; 

4. the initial velocity at which a loop moves after formation; 

5. a, the average length of a new loop after it has formed. 

This last parameter is under very active study at present from both analytical [71, 72, 73] and 
numerical [74, 75, 76, 77] perspectives. It is an important thing to understand because large 
loops live much longer than small loops. Loops lose their length by generating gravitational 
radiation, which is now beginning to be brought under observational limits; since longer 
loops emit gravity waves later, they are more tightly constrained (see ref. [78] for much 
more detail). For the discussion of our numerical results, we first treat them as free and 
independent parameters so as to study how each affects magnetogenesis. On the other hand, 
we take our test values for G/xq from observational constraints: a fiducial value of 2 x 10~^ 
[29, 30], an optimistic value of 7 x 10~^ [31], and the most constrained value of 2 x 10~^ 
[32]; the former two come from combining CMB data with other cosmological observations, 
while the last is the "worst case scenario" limit from pulsar timing. At the end of this 
section, wc present constraints on what we consider to be the best motivated combination of 
parameters, which differs slightly from our fiducial model. Our fiducial model, used wherever 
nothing else is specified, includes: the VOS model for long strings; loop dynamics turned 
on; G/xo = 2 x 10~^; a = 0.01; and vt{t = tp) = 0.1. Where relevant, we have drawn a 
line demarcating the minimum correlation length necessary to seed galactic dynamos and a 
line indicating the minimum magnetic field strength necessary, given a particular dynamic 
amplification time, ^^y- 

4.4 Discussion 

We now discuss the results of the numerical simulations. Unless stated otherwise, all of 
the plots assume the VOS model with loop dynamics, G/xq = 2 x 10~^, a = 0.01, and 
Vt{t = tp) = 0.1. We have systematically explored the variation of our predictions with the 
parameters listed above. 

The scaling of the magnetic field spectrum with G/xq, illustrated in Figures 3, can be 
understood within the context of the calculations in Sections 2 and 3. As G/iq is reduced, each 
loop becomes less effective at generating magnetic fields, and so the magnetic flux is reduced. 
Furthermore, the loop length spectrum peaks at a characteristic length determined by G/xq, 
as given in (3.8) for the OSM model. Since we use only a single scale for loop formation here, 
we find the same result in the VOS model. As G^o is reduced, this characteristic loop size 
falls as well, shifting the magnetic field spectrum to smaller correlation lengths. 



- 19 - 



10-20 

o 
II 

N 

BO 
C 
V 

u 

2 

c 

CO 

10' 10= 10= 10< 10» 10» 

correlation length (pc) at z=10 

Figure 3: Magnetic field strength as a function of the magnetic field's correlation length, for our 
three test values of bare string tension G/io- The plot shows the magnetic field strength at = 10 as 
a function of the proper correlation length at the same redshift. 

The magnetic field strength as a function of correlation length for several values of a 
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Figure 4: Volume of universe suffused with a seed field as a function of seed field flux magnitude. 
The horizontal axis gives the magnetic field strength at z = 10, and the vertical axis shows the fraction 
of the volume of the universe which is suffused with a magnetic field of that flux or greater. 
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Figure 5: Magnetic field strength as a function of the string formation length fraction a. 
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Figure 6: Volume of universe suffused with a seed field as a function of seed field flux magnitude, for 
various values of a. 

is given in Figure 5. The magnetic field spectrum is only weakly dependent on a, thanks 
primarily to two effects. First, while a sets the size of the largest loops, at any fixed time 
the greatest number of loops have characteristic length that is set by GfiQ, as described for 
the OSM model in (3.8). Therefore the peak correlation length and magnetic field flux only 
weakly depend on a. Secondly, loops created in models with larger a will always shrink to 
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Figure 7: Volume of universe suffused with a seed field as a function of seed field flux magnitude, for 
various values of the loop's initial translational velocity, Vt{t = tF)- 

smaller lengths, and thus mimic models in which a is smaller. The peak magnetic flux grows 
slightly with a, since at progressively smaller values of a, the string network must shed energy 
into an ever greater number of ever smaller loops. Thus, the fraction of loops creating fields 
at the peak correlation length grows slowly as a decreases. 

To provide a viable magnetogenesis mechanism, cosmic string loops must not only pro- 
duce sufficiently strong magnetic fields, but must produce them over nearly the entire volume 
of the universe. Figures 4, 6, and 7 explore this issue for different values of G/xq, a, and loop 
vt at formation. In these figures, we have plotted the fraction of a universe suffused with 
a given magnetic field strength or greater. If one chooses a value of magnetic flux required 
to seed a galactic dynamo, the curves give the fraction of the universe in which the string 
network produces the required seed fields. Figure 4 shows the same fall in peak magnetic field 
strength with falling G/xq as does Figure 3. Beyond shifting to smaller values of magnetic 
field strength, the curves look roughly similar: as G/xq is decreased, the loop network contains 
roughly the same number and sizes of loops, which produce the same volume coverage but 
lower magnetic field strengths. Figure 6 shows the same weak dependence of the peak mag- 
netic field strength on a as does Figure 5. The volume coverage falls nearly linearly with a, as 
is to be expected since at smaller values of a the loop network contains fewer large loops. The 
scaling, very roughly speaking, is that the volume swept out by a single loop is proportional 
to a^, while the number of loops produced is proportional to a~^, the combination of which 
gives an approximately linear scaling of volume with a. Finally, in 7, we see an approximately 
linear scaling in volume coverage with loop initial transverse velocity. Faster loops can sweep 
out more volume, but since they remain in any one area for a shorter period of time, they 
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Figure 8: Since loops cliange their size over time, different volumes of space will be endowed with 
seed fields of differing magnitudes. Here we plot a line indicating the percent of the universe covered 
by a magnetic field whose size is greater than or equal to some minimum field strength. Results for 
both network models are shown for comparison, and for each we show results both with (+LD) and 
without loop dynamics. To explain a bit better, an example: the point where the red line crosses 8% 
is at a magnetic field strength of approximately 10~^° G. What this means is that 8% of the volume 
of the universe at the time of galaxy formation was suffused with a seed field whose magnitude was 
equal to, or larger than, 10~^° G. Here we have adopted G/io = 2 x 10~^ and a — 10~^. Only fields 
whose correlation lengths are sufficiently large {Lcorr > 500pc at z = 10) are included. We have also 
included what we believe to be the best motivated model: the VOS model for the long strings, loop 
dynamics, and a = 0.1, but with only 10% of the string network's energy loss going into loops that 
large, with the rest lost to very small loops. 



produce weaker magnetic fields. This is also evident in the plot. 

Figures 8, 9, and 10 explores the effect of loop dynamics and the network model on our 
conclusions. In Figure 8, we plot the horizon volume overswept versus the magnetic field 
strength, as in Figures 4, 6, 7 above. Two effects are evident in these results. The One Scale 
Model covers a slightly large volume of the universe with fields than the Velocity-dependent 
One Scale model, as the OSM tends to overestimate the number of long strings present 
during the matter-radiation transition era, leading to a greater number of loops present 
during magnetogenesis. The inclusion of loop dynamics increases magnetic field strength but 
decreases volume coverage, since loop dynamics decrease loop translational velocity over time. 
In addition to these model comparisons, we also plot the model we believe to have the best 
physical motivation: the VOS model, with loop dynamics included, where string loops are 
formed at 10% of the horizon size, but with only 10% of the loop energy going into loops this 
large. We assume that the remaining 90% of loops are formed near the gravitational radiation 
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Figure 9: Magnetic field strength for the different string models we consider, excluding our best 
motivated physical model (which on this plot is nearly identical with VOS+LD), compared with the 
prediction from long string encounters, where the vortex is created in the region between two oppositely 
moving long strings. 
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Figure 10: For the five model combinations we consider, we plot the fraction of the volume covered 
versus the correlation length. Only fields whose correlation lengths are sufficiently large {Lcorr > 500pc 
at z = 10) are included. 
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back-reaction scale, through loop fragmentation or other effects [78]. In Figure 9, we compare 
the magnetic field strengths in our models with the magnetic fields expected from long string 
encounters. Since the vortices formed by long string encounters are spread over such large 
length scales, the fields generated by those encounters arc concomitantly much weaJier than 
those generated by loops. Finally, in Figure 10, we plot volume overswept versus correlation 
length for these same five models, beginning with the largest correlation lengths. This plot 
includes only field correlations lengths that are sufficiently large {Lcorr > 500pc at z = 10). 

5. Conclusions 

In this paper, we have studied the viability of a cosmic string network for generating the seed 
fields necessary for explaining, after dynamic amplification, the /iG magnetic fields observed 
in spiral galaxies. We have analyzed the behavior of wiggly cosmic strings in the era between 
matter-radiation equality and decoupling. In this era, we find that strings can create vor- 
ticity in the primordial plasma, allowing the Harrison-Rees mechanism for magnetogenesis 
to operate. This effect is, however, at second order in the string's effective Newtonian mass 
density, contrary to earlier estimates in the literature. Though the turbulent eddies in wakes 
behind long strings can form stronger magnetic fields, by our calculations these fields are on 
length scales too small to explain galactic magnetic fields. On the other hand, we find only 
weak support in our calculations for vortices forming in the space between distantly separated 
straight strings, which would have adequate correlation lengths. 

However, by interpolating between these extremes, we have discovered that cosmic string 
loops provide an excellent mechanism for magnetogenesis, with correlation length scales that 
are large enough as well as field strengths much greater than those predicted for the spec- 
ulative inter-string vortices. If we take present work on loop production seriously, then 
approximately 10% of the energy a string network loses goes into strings whose length is ap- 
proximately O.lLff- These loops can oversweep approximately 10% of the universe during the 
relevant epoch between matter-radiation equality and decoupling. With loops preferentially 
oversweeping areas with over densities, this may be just enough coverage. Also, if cosmic 
strings are {p, q) strings from string theory, we can expect a further enhancement of this 
volume coverage by a factor of approximately 3P^ [28]. In those regions that are overswept, 
strings with tensions allowed by present-day cosmological constraints can generate seed mag- 
netic fields large enough to account for today's magnetic fields under the assumption of an 
efficient galactic dynamo for amplification. Although we are making optimistic assumptions 
about these dynamos, we have consistently made conservative estimates for the string-driven 
magnetogenesis. For instance, we are not forced to rely on turbulence or inverse cascades to 
achieve these results; if such mechanisms exist, then the seed fields we predict would only be 
strengthened. 

We achieved these predictions first through careful analytical estimation, then through 
the construction of a thorough computer code. This latter has allowed us to simultaneously 
solve for the evolution of the string network, loop production, loop dynamics, and cosmology. 
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giving us predictions for magnetic fields that take into account a multitude of effects. Beyond 
a bare estimate of the magnetic field amplitude, we correlate the volume of the universe 
overswept by loops with such parameters as G/xq, the string tension, and aLn, the typical 
size of loops. In addition, we compare the results of two popular semi-analytic string network 
models, the original one scale model as well as the more sophisticated velocity-dependent 
one scale model. Because of the versatility supplied by the code, we predict, inter alia, the 
spectrum of magnetic field strength versus correlation length. Though it may be far fetched, if 
future astronomical observations were capable of detecting and characterizing the primordial 
seed field, our predictions are sufficiently concrete that such observations could either confirm 
or refute our model. 

In addition to futuristic observations, our model could also be refuted through the defini- 
tive disproof of efficient galactic dynamos (our model requires < 0.3), or the constraint of 
cosmic string tensions to much less than G/j, ~ 10~®. On the other hand, if cosmic strings are 
observed and possess a tension of around our fiducial value, the case for this magnctogenesis 
mechanism would be greatly strengthened. 

Thus, the presence of a string network might help explaining some large scale magnetic 
fields if efficient galactic dynamos are present. Put in other words, with efficient galactic 
dynamos, we can explain a fj,G with a Gf^. 
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A. Dynamics of rotating loops 

Loops are acted upon by dynamical friction forces as well as the recoil and shrinking from 
the emission of gravitational waves. We must estimate these effects to ensure that the ap- 
proximations we make in Section 2.2 are valid, as well as to properly incorporate them in the 
numerical estimates presented in Section 4.3 and described in detail in Appendix B. In this 
Appendix we derive the differential equations that govern the evolution of the translational 
velocity vt{t), rotational velocity Vr{t) and length £{t) of a loop, and provide approximate 
analytic solutions. 
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A.l Changes in size and shape 

The length of the loop decreases due to gravitational radiation, so that 

m = fralHitp) - TiGfioit - tp) (A.l) 
= 4 - Gr,Mo(i - tp) , (A.2) 

where fr < 1 describes energy loss directly after formation, F; ~ 50 controls the efficiency 
with which the loop emits gravitational radiation [26, 79, 80], io = atp = frOiLHitp), and 
the time scale for loop shrinkage is 

o, 

tshrink = *F = 2000*^ , (-A-3) 

for a = 0.01 and G/iq = 10~^; so a loop with d > 10~^ remains large enough for magnetoge- 
nesis for many Hubble times. 

A further concern is large-scale loop oscillations, which change the shape of the loop 
significantly. (Small scale oscillations of the string are already averaged over to give the 
effective tension T and linear mass density fi). If the oscillation timescale is comparable to 
the rotation period, the rigid loop approximation does not apply. Since the initial velocity 
distribution on the loop has to be quite peculiar to yield a fast oscillation affecting the loop 
as a whole, we expect that these oscillations make only a small fraction of loops unsuitable. 

A.2 Translational Movement 

Three effects determine the drift velocity of a loop: redshifting from Hubble expansion, 
dynamical friction due to dragging of the plasma as computed in (2.18), and recoil from 
gravitational wave emission [81] (often referred to as the rocket effect). The latter causes an 
acceleration of TpGnon/i, where n is a unit vector in the direction of recoil and Tp f« 10 
[81]. (It is fiQ and not \ = ji — T = fi{l — ;Ug//i^) f» 0.56;Uo which determines the emission of 
gravitational waves.) Incorporating these forces leads to 

V, = -H., - + (A.4) 

where t^ is defined in (2.20) and 9min in (2.19). In the matter era, using H = 2/(3i) and 
SirGp = 3H^ gives 

t, = ^, Ci = -GeX. (A.5) 

Since our arguments in the previous subsection show that the loop length £ decreases very 
little over the timescales we are interested in, we take ^ to be a constant. Further, we ingore 
the time dependence in the logarithm and estimate this factor by ^ 

'Svf{tF)tF^ 



« 1^1 )= const, (A.6) 



■''Wc approximate rmax = fvtdt w 3vtt « 3vt{tF)tF, anticipating vt oc 1/a oc initially, due to 

redshifting. 



-27- 



where the time of loop creation is 

tp = e/a. (A.7) 

Different initial values can lead to great differences in the long-term behavior of vt- The 

loop either slows rapidly, or accelerates quickly to relativistic velocities from the rocket effect. 
Below we integrate the dynamical equation (A. 4) in each regime, and derive the limiting 
velocity that separates the two regimes. Without the rocket effect term, equation (A. 4) 
becomes 

St t^v^ ^ ' 

which has the solution 

1/3 



v{{t) = t--l^ 



-3Ci ln(0-Jjt + ^(3Ci ln(^-Jja + vli)^ 



where the last step uses v{tF) = v^^, neglects 3Cialn(6'~^^), and employs the time scale of 
dynamical friction defined by 

t. = (A.IO) 



mm' 



When dynamical friction is irrelevant (A. 4) becomes 

vt = ~vt + C2, (A.ll) 

where C2 = FpGi^o/l and we assume that the recoil is collinear with the velocity. This has 
the solution 

- ^0 f ^V^7l + f ^ V^') , (A.12) 



at J 

where we neglect MC2 in the last step and introduce the relevant time scale tr for the rocket 
effect 

Setting tr = tf and solving for the initial velocity yields 

/x5\ 1/12 X . N 1/4 

vum = G^/'[^) (^) 2V-a^/^(ln^-L)^/-, (A.14) 
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Hence, for vt -C vnm dynamical friction is more important, whereas for larger velocities, the 
rocket effect predominates in the long run. 

What is a reasonable initial value for vt! If large loops are created in the matter era, we 
expect their velocities to be comparable to the RMS velocity in the string network vrms ~ 0.6 
^. If we average over the small scale wiggles, the RMS velocity drops down to vrms ~ 0.15 
[41], so the initial translational velocity of the loops should be smaller too. Rotating loops 
share kinetic energy between translation and rotation, so an initial value of vt ^ Vr, say down 
to 0.4, may be reasonable. This velocity is still larger than the limit velocity (A. 14), even for 
the largest loops produced in the network: for a = 0.01 and the bound GA/0.56 = G/xq ^ 
10"''', the limiting velocity becomes < 0.0072, where we approximated ln0~^^ ~ 19, 
based on (A. 6) with vt{tp) = 0.4. Hence, the translational velocity of loops which are created 
in the matter era is given by 

for all relevant values of a. The loop initially slows due to redshifting, even though the rocket 
effect dominates over dynamical friction. Indeed, for a ^ 0.01, vq = 0.4 and TfiQ = 10~^ we 
have tr ~ 197£/q; = 197 ip. In this case, vt drops down to the minimal value of vt ~ 0.024 
at t pa tr3/4, and it increases linearly thereafter up until loop shrinking becomes important. 
Since it takes more than 1000 for a loop to accelerate so that it moves faster than vq again, 
we expect the majority of loops to have a translational velocities of order vt ~ 0(1O~^), which 
we use as a rough estimate of vt for loops created between teq and tdec- Since tr < tghrink 
from (A. 3), our assumption of a fixed loop length £ is justified. 

A. 3 Rotational Movement 

The rotational velocity of a loop is influenced by three effects: dynamical friction from 
plasma drag, emission of gravitational radiation ^ which produces a torque [82] 

Tgr = -iG/llVgr (A. 16) 

where Tgr ~ 5 [82], and the shrinking of the loop from gravitational wave energy emission. 
The torque Tdrag from dynamical friction is computed using the arguments of Section 2.2, 
applied to a single loop rotation. During one rotation, spanning a time At ~ £/vr, the 
surrounding plasma acquires an angular momentum of roughly 



A J ~ lAppiasma ~ ^^pVx ~ (A. 17) 



^Recent simulations and analytic arguments suggest two distinct classes of loops [72, 78]: small, highly 
relativistic loops, and large ones, which have the velocities we consider. 

^There is no rocket effect for angular momentum: numerical studies show that the emission of gravitational 
waves always decreases the angular momentum [82], oven though a rigorous mathematical proof is lacking. 
The fundamental mass density no determines this effect and not the effective mass density A. 
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which gives the torque 



Tdrag ~ " V —^ P' (^'l^) 



7 

Loop shrinkage enters via the expression for the total angular momentum J, 



J 



j = A (2eivr + fvr) , (A. 19) 

with £{t) from (A. 2). We neglect redshifting of the rotational velocity for loops well within 
the horizon. 

Unlike the translational velocity, it is possible to obtain analytic solutions for the rota- 
tional velocity with dynamical i{t). Dynamical friction is negligible, since the ratio of torques 
is 

= ^ o „ ^ « 1-3 X 10"^^ < 1 (A.20) 



'gr 



7 nlv^VgrGTT ■ 



where we use SirGp = 4/(3t^) as well as X/po ~ 0.6, a ^ 0.01, Vr ~ 0.4 and Fgr ~ 5. This 
leaves J = Tgr, leading to 

_ 2GTipoVr{t)-Cgr 2n 

'''•-4-Gr,Mo(t-M' ^ ' 



where 



Cgr = ^^^^ , (A.22) 



and i{t) is given by (A. 2). Assuming the initial condition Vritp) = vq, then Vr{t) is 

1 Cgrfjt) + 2voa 

2 GFipofit) + a 



^.(*) = i%^%^ (A.23) 



with 



2 



/(t) = 2«(l- Aj +Gr;Mo(l- Aj . (A.24) 

As with the translational velocity Vt, there is a critical initial rotational velocity that 
determines the future evolution of Vr, given by vj-^^ = Cgr/ {2GTiPq). For the parameters in 
table 2 we get > 1, indicating that Vr decreases for all loops under consideration. 

What is the initial velocity of a loop? Loops can be produced by self-intersections of 
a single string. In this case, we expect the intersecting pieces to move roughly in the same 
direction. As a consequence, the majority of the energy will go into translational movement 
and a loop with very little angular momentum but large momentum results, so that Vj- <^ Vf. 
However, Loops can also be created when two strings, moving in opposite directions, intersect 
and chop off a loop. Loops formed in this way have very little momentum, but large anglar 
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momentum, so that initial rotational velocities comparable to the RMS velocity of the network 
result Vr vrmSj while Vf <^ Vr- 

An ideal loop for magnetogenesis lies somewhere in between. Optimally efficient loop 
magnetogenesis requires rotation to stir up the plasma, but also rapid translational motion 
to oversweep a large fraction of the universe. To account for kinetic energy in the form of 
the small scale wiggles and oscillations of the loop, as well as the angular momentum that is 
radiated away in gravitational waves immediately after its creation, we take Vr ~ 0.4 < vrms 
to be a conservative initial value. 

B. Code Implementation 

To fully account for the background cosmological model, the forces and torques on rotating 
loops, and the behavior of the string/loop network when the universe is not in a scaling 
regime, we employ a computer program which implements the full set of differential equations 
described in this paper. Here we discuss the operation of this code. 

The code follows an array of loop cohorts over time, while simultaneously tracking some 
averaged quantites describing the cosmological model and network of long strings. At each 
time step, the code updates the background cosmological quantities (Hubble parameter, en- 
ergy density components, etc.) according to the Priedmann and energy conservation equa- 
tions. It also updates the properties of the long string network by using either the one-scale 
model (Section 3.1) or the velocity-dependent one-scale model (Section 3.2). Both the OSM 
and VOS models include a production rate for loops at a size determined by the background 
cosmological parameters, so in each time step the appropriate number of loops are added to 
the corresponding cohort in the array. 

In each time step we include all of the physics of the loop dynamics and magnetic field 
generation for each loop cohort. Each individual cohort collects all the loops created during a 
block of time steps. We typically use 10^ steps equally spaced in log(t) from rcdshifts z ~ 10^ 
until the present, with blocks consisting of 10^ steps each. The loops created during each time 
block have similar properties, but there are slight differences between the loops added at the 
beginning and end of a single block. The code accounts for these differences by updating some 
of the cohort parameters through a weighted average of the properties of the newly added 
loops and the properties of the loops already in the cohort. Each cohort tracks a number of 
different variables: 

1. The comoving number density of loops in the cohort, as computed by the OSM (3.2) 
or VOS equations (3.12). The comoving number density is not averaged at each time 
step, but accumulates new loops as they are created. 

2. The physical length of the loops in the cohort, which decreases due to the emission of 
gravitational radiation, with the rate given by (3.4). This is averaged over the cohort. 

3. The translational and rotational velocity of the loops, which are computed using (A. 4) 
and (A. 19), respectively. This includes Hubble damping, the rocket effect, dynamical 
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friction, torque from graviatational wave emission, and the shrinking moment of inertia 
of the loops. These properties are averaged over the cohort. 

4. The fractional comoving volume overswept by the loops in the cohort, which depends on 
the comoving number density and translational velocity of the loops. This is a proxy for 
the fraction of the universe over which the loops in the cohort may generate magnetic 
fields, and is accumulated (not averaged) as new loops are added. 

At each time step, the program updates each of these quantities, from the time the first 
loops are added to the cohort until the cohort has completely evaporated into gravitational 
radiation. 

The ultimate objective is to calculate the magnetic field strength as a function of scale 
created by the loops. This is tabulated in an array which tracks the comoving magnetic fluxes 
a^{t)B{t) on a selection of length scales. Each magnetic field bin also tracks the fractional 
comoving volume of the universe which contains magnetic fields at that length scale. An 
individual loop cohort generates magnetic fields on large comoving scales when the loops are 
newly created, and then on progressively smaller scales as the loops shrink. Furthermore, the 
magnetic field at each comoving correlation length gets contributions from many different loop 
cohorts over time. Therefore the array tracking magnetic field strengths is indexed differently 
from the one tracking the loop cohorts, with each array element corresponding to a fixed 
range of comoving lengths. 

At each time step, the code steps through the loop cohorts, and assigns the magnetic fields 
generated by each cohort in that time step to a bin in the magnetic field array, determined by 
the size of the loops in the cohort. The comoving magnetic flux is then updated by averaging 
the flux already present and the flux created by the loops in the cohort, weighted by the 
comoving volume traversed by the loops in the cohort and the total comoving volume already 
present in the bin. The comoving volume represented by the bin is incremented by the volume 
traversed by the loop cohort over the last time step. At the end of the simulation each bin in 
the magnetic field array represents the average magnetic field strength in a certain fraction 
of the universe over which the field is correlated over a given length. 

To compare the field created by the loop population to that created by the long strings, 
a parallel calculation tracks the magnetic field created by the long string population. This 
is much simpler than the loop calculation, primarily because the correlation length for long- 
string magnetic fields is given by the interstring separation, which is usually some nearly 
constant fraction of the horizon length. Thus each comoving correlation length only gets 
contributions from a fixed moment in cosmic history. 
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Parameter Value 



Source 



description, first used in (equation) 



a 

fr 

a 



gr 



r 



dy 



^dec 
^9f 

ieq 
^dec 

to 

vrms 
vrms 

h 

P 
P 



ci 

C3 



2 X 10-7 

1.9/xo 

1.5//0 

0.01 

0.7 

fraH-\tF)/tF 

50 (★) 

50 < < 100 

45 < r, < 55 

50, 80 

5 

10 

0.2 Gyr (*) 
0.2 < r^jJ/Gyr < 0.8 
r^;>(l-l-l-4)Gyr 



dy 
dy 



= 0.3 Gyr 

r^^i = 0.5Gyr 
r-; = 2.2Gyr 
1089 
6 

10 w 

1.6 X 10^2 g ^ 51 kyj. 

1.2 X 10^3 g ^ 3gQ kyj. 

13.7 Gyr 

0.60 

0.15 

0.1277 

0.732 

1/2 < /3 < 1 
1 (*) 

10-1 < P < 1 
10-3 < P < 1 
0.21 (0.2475) 
0.18 (0.3675) 
0.28 



[29, 30] 

[45], radiation era 
[45], matter era 
[26, 45] 

[87, 26] 

Defined in (4.13) 

[79] 
[80] 

[26, 82, 84] 
[82] 

[83, 81, 82] 

[86, 60] 

[15] 

[68] 

[64] 

[67] 

[85] 

[22] 



[85] 

[41], matter era 

[41], matter era 

[85] 

[85] 

[23] 

cosmic strings 
D-Strings [53] 
F-Strings [53] 
[27, 28] 



string mass/length 
effective mass/length (2.3) 

size of large loop/horizon (3.4) 

loop redshifting energy loss (3.4) 
new loop length over formation time 
efficiency of grav. wave emission [43] (3.4) 



grav. wave emission yielding torque (A. 16) 

rocket effect (A. 4) 
dynamo efficiency (4.5) 



redshift at decoupling 
redshift of galaxy formation 

time of matter-radiation equality 

time of decoupling (4.10) 

age of the universe 

RMS velocity in network 

RMS velocity avg. over corr. length 

matter fraction (4.7) 

Hubble parameter (4.7) 

(2.7) 

intercommutation probability 



VOS parameters in radiation (matter) era 



Table 2: Our parameters. When several values are given, we select the one marked by a (i*r). 
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